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Abstract. A bordering of GUE matrices is considered, in which the bordered row consists 
of zero mean complex Gaussians N[0, cr/2] + jN[0, cr/2] off the diagonal, and the real Gaussian 
N[/i, a/\/2] on the diagonal. We compute the explicit form of the eigenvalue probability 
function for such matrices, as well as that for matrices obtained by repeating the bordering. 
The correlations are in general determinantal, and in the single bordering case the explicit 
form of the correlation kernel is computed. In the large N limit it is shown that fi and/or 
a can be tuned to induce a separation of the largest eigenvalue. This effect is shown to be 
CNj ■ controlled by a single parameter, universal correlation kernel. 
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As is well known, in what was the hrst application of random matrix theory to physics, Wigner 
introduced a particular ensemble of real symmetric matrices as a model of the highly excited 
states of complex nuclei (see e.g. [H]). Thus the features of the Hamiltonian hypothesized to be 
responsible for the distribution of these states — in a minimalist reduction chosen to be simply 
f^ time reversal symmetry and the absence of a preferential basis — were observed to be features 

shared by the ensemble of real symmetric matrices with probability density function (PDF) 
proportional to the Gaussian exp(— TrX 2 /2). Hermitian matrices of size N x N, with N large 
(formally N — > oo) are chosen because the discrete portion of the spectrum consisting of an 
infinite number of levels are being modelled. The presence of a time reversal symmetry means 
the basis elements of the Hamiltonian and thus the Hermitian matrix can be chosen to be real, 
iq ■ thus explaining the use of real symmetric matrices, and the PDF is the simplest that is invariant 

under mappings X i— > RXR where R is a real orthogonal matrix, this having the physical 
interpretation of there being no preferential basis. This invariance lends the name Gaussian 
orthogonal ensemble (GOE) to this class of random real symmetric matrices. 

Of direct relevance to the concerns of the present paper is a variant of the GOE specified 
by the PDF proportional to exp(— Tr (X — Xo) 2 /2). Here Xo is a fixed real symmetric matrix 
specifying the mean of the Gaussian distribution. Scaling this PDF by a parameter t so that it 
reads exp(— Tr (X — Xo) 2 /2£), one obtains a well known (see e.g. [T31ITT]) realization of the Dyson 
Brownian motion model [7] . In this the corresponding eigenvalue PDF evolves in t according to a 
certain Fokker-Planck equation, with the initial condition that the eigenvalue PDF is a product 
5_j ■ of delta functions at the eigenvalues of Xo. 

Shifting the mean in the PDF for the GOE is equivalent to forming an ensemble of matrices 
of the form X + Xq where X is a member of the GOE and Xq is as above. Special to this class 
of matrices, and a primary concern of our work herein, is the case that Xq is of low rank. The 
simplest of these is when Xo = ee\e[ , where ei is the elementary column vector with 1 in the first 
entry. By orthogonal invariance of the distribution of X the corresponding ensemble of matrices 
has identical spectral properties to the shifted mean GOE ensemble in which Xo has all entries 
equal to e. This ensemble appeared in an analysis of a spherical spin glass due to Kosterlitz, 
Thouless and Jones [T3] (see also pQ), and it has also been used as a model Hamiltonian in the 
study of mesoscopic quantum structures 19]. 
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Tuning the parameter e by writing e — c/y2N, (c > 0), it was shown in |14j that in in the 
large N limit the largest eigenvalue undergoes a phase transition as a function of c. Thus for 
< c < 1 the leading order location of the largest eigenvalue is at y/2N (which is the edge of 
the support of the spectrum for e = 0), while for c > 1 it is at 

(1.1) (N/2)^(c+l/c) 

and thus separates from the leading support of the spectrum. The same effect occurs for the 
ensemble of complex Hermitian Gaussian matrices with PDF proportional to exp(— TrX 2 ). This 
specifies the Gaussian unitary ensemble (GUE). Moreover, is this case the analytic form of the 
scaled correlation functions can be computed exactly jTTJ |S] . 

In this work we will analyze the scaling regime of this eigenvalue separation effect for a 
generalization of the additive perturbation X + Xq, X a member of the GUE and Xq of low 
rank. The main generalization to be considered is the class of bordered matrices 



(1.2) 



c N v N 
v N G N 
d 



where Gm a member of the GUE, cm = N[//, a/\/2\ and vn a N x 1 vector of complex Gaussians 
N[0, cr/2] + iN[0,ct/2]. In the case that a = 1 this corresponds to perturbed GUE matrices 
X + diag [fi, r ] with X a member of the (N + 1) x (N + 1) GUE and the notation 0™ denoting 
the eigenvalue repeated n times. Up to a similarity transformation, such matrices are equivalent 
to (N + 1) x (N+ 1) shifted mean GUE matrices X + ee\e{, e = fi/(N + l). We will furthermore 
consider interations of (jl.2p obtained by bordering it with a new first row and column with 
elements distributed as in the previous first row and column. Thus we are generalizing the 
shifted mean GUE ensemble by also allowing for a different variance along the first row and first 
column, or more generally the first r rows and columns. 

In a previous paper [3] we gave an extended discussion of the relevance of matrix ensembles 
with varying mean and variance to contemporary studies in mathematical statistics and applied 
mathematics. Most notably these include the study of spiked models in the analysis of multivari- 
ate data [2], and stability questions relating to ecological webs [15 . Some universal behaviour 
has been found. Thus it has been shown [17l |6] that the critical region of the r — 1 complex 
spiked model (complex Wishart matrices with a covariance matrix having a single eigenvalue 
different from 1) is the same as that for the rank 1 shifted mean GUE. We will exhibit the same 
critical correlations for the model (jl.2p as a function of a in the case /i = 0, or more generally as 
a function of both \x and <r, when tuned about the eigenvalue separation point. 

Our first task, addressed in Section 2, is to give a determinantal formula for the eigenvalue 
PDF of matrices (jl.2p and their extensions defined by bordering further rows and columns. The 
determinantal structure is used in Section 3 to compute the correlation functions for fixed values 
of the parameters. With the parameters tuned in the neighbourhood of the eigenvalue separation 
point, the N —¥ oo scaling limit is studied in Section 4. 

2. Joint eigenvalue PDF 

We begin by observing that the eigenvalue PDF for matrices of the form (|1.2[) is unchanged 
by the replacement of Gn by diag Gat. 

Lemma 2.1. Consider matrices of the form Al.ty) . generalized so that Gn is any Hermitian 
matrix independent of cjv and vm ■ Then il.ty) and the random matrices 



(2-1) 

have the same eigenvalue PDF. 



CN 

v N diag G N 
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Proof. Since Gn is assumed Hermitian we can write U'GnU = diagGjv for some N x N unitary 
matrix U. Noting too that a complex Gaussian is unchanged in distribution when acted on by a 
unitary matrix we see that (I2.1J) has the same distribution as 



1 



TV 



CN % 

v N G N 



1 



A' 



U 



□ 



Oat W 
This implies the result. 
Suppose 

(2.2) diagG/v = diag[ai,a 2 , . . . ,a w ] 

and let us denote by x(A ^~ a ) the indicator function for the interlaced region 

(2.3) Ai > oi > A2 > &2 > ■ ■ ■ on > Ajv+i. 

We can use theory from [T^] to give the explicit form of the eigenvalue PDF of l|2.1|) for given 
{ai}. 

Proposition 1. Consider random matrices 1)2.1)) . with diagGjy given by 12.2(1 and cm , Vn having 
Gaussian distribution specified below 111.2)) . With {a^} given, the eigenvalue PDF is 

/ V / 2\ 2W + 2 n~ 



V a ) 



-2 v-^ iv + 1 \ 2 



E j= i A^-Ci-o- )E j= i« Je - 



M"Y 



(2.4) 



xe 



a/wr-'fE^S'^-E^i-j 



.m 



<i<fe<jv+ 



1 (\ j - A fe ) 



a-k 



x(A>-o). 



lll<j<fc<w( a j 

Proof. Elementary manipulation of the corresponding characteristic polynomial corresponding 
to (|2.1[) shows that the condition for A to be an eigenvalue is 



(2.5) 



= A — cn 



N 



4 ] \ 2 



where v N denotes the jth component of vn ■ By the specification of the latter given below (|1.2|) 
we see that 

\v$\*±r[l,o*/2]±(o*/2)T[l,l], 

where T[k, c] refers to the gamma distribution, specified by the PDF proportional to (s/c) k ~ l e~ s /° 
Writing 



then allows (J2.5I) to be written 


C N = {-7=)C, Clj 


(2.6) = 


N 
= \ c V."* 

. A — a, 



where lOj =T[1, 1]. 

The distribution of the zeros of the zeros of the random rational function (|2.6p . and thus 
the eigenvalues {\j}j=i,...,n hi the case of cjv fixed, is given by [ITj (Prop. 8 with n = N 7 
si = ■ ■ ■ = sjv = 1, A n. A — c, aj i— >• 5j — c). Explicitly, in terms of the original variables {%■}, 
{\j} we have that the PDF for the zeros is supported on the region x(A >- a) and is subject to 
the further constraint 

JV+l TV 

(2.7) ^2\i = Y,ai + c N . 



1=1 



2 .2 



4 K.E. BASSLER, P.J. FORRESTER, AND N.E. FRANKEL 

In this region the PDF has the explicit functional form 

( 2 g) (V2\ 2N+1 Tll< 3 <k<N + l(^ - A ^ -.-^S'a;-^, a*) e a- 

^ ° ' Ul< j<k <N( a ] - a k) 

In the case that c = N[/x, <r/\/2], the reasoning of [HI Corollary 4] tells us that the constraint 
(|2.7p should be removed, and in (12.8[) the term exp(— a~ 2 c N ) should be replaced by 



(2.9) ,/4^- 2 ES 1 *,-E£ 1 ,) e --V ! . 



7TCr 2 



This gives (g3). D 

Consider now a random Hermitian matrix Gn+t specified so that the first r diagonal entries 
have distribution N[/x, a/\^2], the entries to the right of these rows have distribution N[0, a/2] + 
iN[0,cr/2], while the bottom N x N sub-block is a CUE matrix. According to Lemma [2.11 we 
have the inductive construction 



(2.10) G 



N+s 



CN+s-1 %"+ s -l 
VN+s-1 Gn+s-1 



(s = l,...,r). 



We can use Proposition [T] to give the joint PDF of the eigenvalues {A'™ +2 '} s= o ! i...., r , A^ :— 
(Ai, . . . , Ajv). 

Proposition 2. Let Gm be a member of the GUE, and let Gn+s (s = 1, . . . , r) be specified by 
112.1 0\) . The joint eigenvalue PDF of these matrices is given by 

1 (20- 2 ) rN+r2 / 2 ^ a -2 E « + - (A (« + ") ) 2 +2M(T -2 j^+, a (» + .) 



-e 



C N (27^/2 

n (Af +r) - \{ N+r) ) n x(a (w+s) y a( jv+s - i > 



X 

l<j<k<N+r 8=1 

(2.ii) x [] (xW-xW^+'-'^uirr-^-'^-iV 

l<j<k<N 



where 



N 

I 



(2.12) C N = 7r Ar / 2 2- Ar ( Ar - 1 )/ 2 []. ? 

3=0 

With {Qj{x)}j=o.i,..., {Pj{ x )}j=o,i,... two sets of arbitrary monic polynomials labelled by their 
degree, this can be rewritten to read 

1 (2a- 2 Y N+r2 / 2 -2 E « r(A («+") )2+2 ^-2 E « r A <« + ,> 
C,v (2ttV/2 

r 

x det[g i _i(A{ c JV+r) ] j ^ = i i ... i Ar +r JJ det[x 



<N + s ) >A (N + s -l)Jj ! fc = l,...,Af + S 
S — 1 



(2.13) xdetbj _ l(A W kfc= ,^ Are (-i + - 2 )E^ 1 (Ar , ) 2 —- 2 E^ 1 ^ 



Wftere <7 W = (-l)N(N-l)/2+(N+r)(N+r~l)/2 C 



N- 



Proof. The conditional eigenvalue PDF of Gn+s, given the eigenvalues of Gn+ s -i is equal to 
(|2.4p with A = \( N+S \ a = \( N+S ~ 1 \ Forming the product over s (s = r, r — 1, . . . , 1) gives the 
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joint eigenvalue PDF of {G\zv+s}s=i,...,r given the eigenvalues of Gm- But the PDF for the latter 
is equal to [TTJ Prop. 1.3.4 multiplied by AM to account for the ordering of X^] 

(2-14) ^f[e-^ 2 II (A^-Af ) 2 , 

N j=l l<j<k<N 

which we multiply this conditional probability by to get (|2.11[) . The form (|2.13l) follows from 
(|2.11j) by noting that 

(2-15) x(A y a) = det[xA J >aJj,fe=i,...,w+i 

valid for ai > a-i > ■■■ > cin > aw+i := — oo (see e.g. [TTJ Prop. 5.9.1]), and that, with 
{ r jW}j=o,...,N a set of monic polynomials, 

(2.16) [] ( a i - fl fe) = (-l) N{N - 1)/2 Mrj-Mhk=i,...,N. 

l<j<k<N 

which is a consequence of the Vandermonde identity. □ 

Our present interest is in the distribution of \ t - N+r '> only. Thus we must integrate (|2.13[) over 
A (w+ S ) ; s = o,l,...,r-l. 

Proposition 3. Let G^+r be defined by the recursive construction H2.10)) , and let A' Ar+r ' denote 
the N + r ordered eigenvalues. The PDF of \( NJrr > is equal to 



(2.17) x detfe_ 1 (Ai Ar+r) ] J , fc=1 ,..., JV+r det \[h k -i, r -i(xf +r) )]^_,...^ + r [(\f +r) ) r - s }^_ 



N + r 
k = l.....N uv J J s=l.....r 



where 

(2.18) h k , r {x) = -. I (x-u) r e(- 1+a ^ u2 - 2 ^~ 2u Pk (u)du. 



I 

' • J — oo 

Proof. We integrate over \^ N \ \^ N+1 \ . . . , \( N + r ~ l ) in order. First note that the integrand is 
symmetric in A*-^ so the corresponding integration can be taken over all M. N provided we divide 
by AH. It is also true that both determinants involving A^' are individually anti-symmetric in 
A^. One of these is det[pk-i(Xj ]j,k=i,...,N- The latter property means we can replace this by 
its diagonal term A ? "!T^[ fc _ 1 Pk-i(XS ). The A^ dependent terms in the integrand are thus a 
product over factors involving X k , times the other determinant. The former can therefore be 

multiplied into column k (k = 1, . . . , N) of the latter and the integration over each A^. can be 
done column-by-column to give 

(2.19) det[[Vo(A^ +1) )],=i n + i [1} 3 =i,...,n+i 

k=l,. . . , N 

Here we have used the fact that X A («+i) >A (N) = 1 for k = N + 1 since \ N ^ 
The quantity (J2.19I) must be multiplied by 



det 



Xa( jv+2) >a< 1v + 1) 



j,fc=l,...,JV+2 



and the variables A( w+1 ) integrated over. We make use of an analogous strategy to that just 
used. Thus we note that the integration over \( N+i > can be taken over all R^" 1 " 1 provided we 
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divide by (N + 1)! and furthermore (J2.19I) can be replaced by 

(N + iy.f[h k Axi N+1) ). 



fc=i 



Multiplying this into the columns of the remaining a:^ Ar+1 ^-dependent determinant, and integrat- 
ing over \( N+1 > column-by-column gives 



(2.20) 



det 



h k ,o(v)dv 



j = l,...,JV + 2 
fe=1.....2V 



5,- Jj - =l,...,JV+2 [ i -\j=l,...,N+2 



Here, to obtain the second last column we have used J du = x + c and have subtracted c by using 
the fact that the final column is all l's. Furthermore, the integral in (|2.20p can be simplified 



h kfl (v)dv 



-r-{v — x)hk,o(v) dv = hk,i(x), 

av 



where the second equality follows by integration by parts. 

Continuing this procedure until we have integrated over all A^, \( N+1 \ . . . , A^ JV+r_1 ^ gives 
([2~TTj) for the marginal distribution of A (JV+r) . D 

Some checks on (|2.17|) are possible. Thus the case N — corresponds to an r x r GUE type 
matrix with joint distribution of elements proportional to exp(— cr~ 2 Tr X 2 ). We find that indeed 
(|2.17l) with N = reclaims (12.14)) with N i-> r, x^ i-> xj /a as required. Furthermore, with 
r = we should reclaim what we started with — the N x N GUE. Using the general formula 

1 



lim 
r-»o+ (r — 1) 



(x-u) r - L f(u)du = f(x) 



we can check that (|2.17[) so simplifies. 

The case a = 1 corresponds to an additive finite rank perturbation of an (N+r) x (N + r) GUE 
matrix. Thus extending the remark below (jl.2p we then have Gm+t = X + /idiag [l^O^] The 
eigenvalue PDF in this case can be evaluated using the well known Harish-Chandra/Itzykson- 
Zuber integral [TTJ |B]. First, to make sense of (|2.17[) in the case a = 1, we must subtract a 
multiple of the final column from the first N columns to replace the lower terminal of integration 
therein by (otherwise the integral is not convergent). Then setting 



(2-21) ^_!(^) = ( - — ) e 

and integrating by parts shows 

[r- !) ! Jo 



1 \ r rir 



du r 



1 \ r 



2// 



(e 



-2/ia; fc-l 



Sk,l)- 



It follows that, up to proportionality, (J2.17I) simplifies to 



[{Aj ) \j=l s ..,N+r 



x det 
in keeping with the known expression. 



N+r 

Mr r \xf^y 



j = l,..,W+r 
s = l,..-,r 
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3. Correlations for r = 1 

3.1. The correlation kernel. The simplest case of (|2.17[) beyond the GUE itself is r = 1. Here 
we take up the task of calculating the corresponding correlation functions. Our strategy is to 
seek a choice of {qj(X}}, {pj{u)} such that the structure of a biorthogonal ensemble [16j[5], jTTJ 
Section 5.8] 

1 N 

(3.1) -^— T\ W2(\i) det[£j(\k)]j,k=i N+idet[r)j(\ k )]j,k=i,...,N+i- 

Cn f=i 

Once this is achieved, the general fc-point correlation is given in terms of a k x k determinant 
according to 

(3.2) /9 ( fc)(Ai,...,Afc) = det[K N+ i(\j,\i)] Jt i=i t ... !k 
where the so called correlation kernel K^+i is given by 

JV+l 

(3.3) K N+1 {\,v) := (w 2 (\)w 2 { V )) 1 ' 2 £ (g^+ik.a&M'fcM- 

p,q=l 

In (|3.3[) g^_, i refers to the inverse of the (N + 1) x (N + 1) matrix of inner products 

(3.4) g N+1 = / w 2 (X)^(X)r lk (X)dX 

LJ-oo Ji,fc=l,...,JV+l 

We will seek a form (|3.ip so that the inverse of Q3.4[l is simple to compute. 

Proposition 4. Lei {_ff p (a;)} denote the Hermite polynomials, characterized by the orthogonality 

/oo 
e- x2 H j (x)H k (x) dx = TV^j- fe , A/} = 2 J j!Vr?. 

-OO 

Introduce {ci p }, {/3 p } as expansion coefficients 

: / e (-l+a-^-2^--u du = y apffp(l) 
J p=0 

OO 

(3.6) e {i-„-*)x>+2^ x = Y^PpH^x). 

p=0 

Write 

N-l 

(3.7) e "(*W)/2 J- J—H^^H^iy) =: K^(x,y) 
(this is the correlation kernel for the (N — 1) x (TV — 1) GUE) and set 

( 3 - 8 ) Vj(x) = < oo . 

The correlation kernel li3.S\) for the r = 1 case of \2.1Tty is given by 

K N+l (x,y) = K%™(x,y) + "' - " f (J™H N _ x (x) - ^H N (x)) m (y) 

(3.9) + (^H N (x) - J^-H N ^(x))r, N+1 (y)\ 
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Proof. For notational convenience write 

z{u) = e (-l+^ 2 )« 2 -2M^ 2 «_ 

In terms of this notation, specify the monk polynomials {pj(x)} in (J2.17I) according to 

(3.10) ^)--jYB^)i^i^ x)H ^ x) ) c*=i.-,tf). *(*) 

(cf. (|2.2ip ). Furthermore, specify the monic polynomials {qj(x)} by 

q j (x) = 2-iH j (x) (j=0,...,iV). 
Noting that the first of the specifications p. 101) implies 

2"J 



1 



z(u)pk(u) du 



—z^H^iXj) (k = l,...,N) 



we see that (12.171) with the above simplifications and r = 1 is proportional to 



1 _ sr^N + l ,2 

e ^=i ^detli/fe-iCAj)]^^! 7v+r 

Oat 



(3.11) x det 



*(A,0 J- 



z(u) du 



j=l,...,JV+l 



[-fffe-i(Aj)] j=i,....jv+i 



fc=l.....JV — 1 



*(*i) 



j=l,...,JV+l 



where, with Cm is specified by (12.121) . 

Cn = Otv2 <j \cr — l) 7t ' 



and we have set A,- ' = A,-. 

The expression (|3.11[) is the form (|3.ip of a biorthogonal ensemble with 

&(*) = H j-i( x ) 

f H^x), j = l,...,N-l 

{ E£Ljv-i a P H vi x )i j = N + l 

and gAr+i a diagonal matrix for rows j = 1, . . . , N — 1 with entries A/}, and a 2 x 2 block diagonal 
matrix with respect to the last two rows and columns, with 2x2 block 



0N-1 OiN-l 

(3 N a N 



(3.12) 

It is thus straightforward to calculate gjy+i, which when substituted in (13.3[) gives (J3.9D . D 



3.2. The expansion coefficients. We seek the evaluation of the expansion coefficients in (|3.6[) . 
Since 



(3.13) Afj 



putp 



dxe~ a x +2 ^ x H p (x) / due(- 1+ ° > ~ 2 ^ 



OO 



,-2\„,2 n,,„—'2„ 



/OO 
e-^ 2x2+2 ^ 2x H p (x)dx 
-OO 

this requires us computing some integrals. 
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Proposition 5. We have 

(3.15) M p f3 p = e^ 2 (l - a 2 f/ 2 tf p ( M (l - a 2 )" 1 / 2 ), 

(3.16) Af p a p = d(l - ^/^(Ml - a 2 )- 1 ^) + C2 (l _ CT 2 f/ 2 /i p ( M (l - a 2 )" 1 / 2 ), 
lyftere 



h p (x):= / e - gp( " } du 



and 



..-Ma. Mtfi-0-"V 



(1 - a 2 )V2)/ ll ( Ai (i _ CT 2)i/2) _ 2ixh (n(l - a 2 )- 1 / 2 ) ' 

= a 2 

C2 ~ (1 - a 2 ) 1 / 2 )/i 1 (/u(l - a 3 )- 1 /*) - 2nh (n(l - a 2 )- 1 / 2 ) ' 
In the case \i = i/iese simplify to give 

(3.17) A/" 2p a 2p = 0, A/^+i&p+i = 0, 

(3.18) A/- 2p+1 a 2p+1 - v^2 2p p!aV 2 - 1) P , 

(3.19) AA 2p /3 2p = ^<r(. 2 - 1)". 

p! 

Proof. The evaluation (13. 15[) follows immediately from the well known, and readily verified, 
integral identity 

e -^-vf/^H n {x) dx = (2ttu) 1/2 (1 - 2u) n/2 H n (y(l - 2u)~ 1/2 ). 

Alternatively, we can make use of the three term recurrence 

(3.20) H p+1 (x) = 2xH p (x)-2pH p _ l (x) (p = 0,l,...) 
and the differentiation formula 

(3.21) -^-H p (x) = 2pH p _ 1 (x) 
to deduce the recurrence 

(3.22) AA p /3 p = 2/iAA p _ 1 /3 p _ 1 +2(p-l)(a 2 -l)AA p _ 2 /3 p _ 2 0=1,2,...). 
With N P P P = (1 - o- 2 y/ 2 y p this reads 

(3.23) y p = 2/i(l - iT 2 )- 1 / 2 ^-! - 2(p - l)y p _ 2 . 

Now ({3T2"3l with p >->• p + 1 and x = p(l - cr 2 )" 1 / 2 is identical to (J3T2U| . This establishes (|3~15)) 
up to the value of Ao/3o, which can be checked directly. 

The advantage in setting up a recurrence is that the same strategy works for the integral 
(13.141) . Thus making use of (J3.20I) and (|3.21[) and integration by parts we deduce 

(3.24) N p a p = 2^A/" p _ia p _i + 2(p - l)(cr 2 - l)A/" p _ 2 a p _ 2 + er 2 V^ P ,i (p = 1, 2, . . .), 
This with M p a p = (1 — a 2 ) p ^ 2 x p reads 

x p = 2/*(l - a 2 )- 1/2 x p _i - 2(p - l)x p _ 2 + ct 2 (1 - o- 2 r 1/2 <5 p4 . 

Hence {£ P } P=2 ,3,... satisfies the recurrence (|3.20|) . with x = fx(l — a 2 )^ 1 / 2 and for p = 1,2, 

The recurrence, being of second order, has two linearly independent solutions H p {x) and h p {x) 
so we have 

(3.25) A/> p = ci(l - a 2 y' 2 H p (n{l - a 2 )- 1 ' 2 ) + c 2 (l - a 2 )*'/ 2 /^! - <rT 1/2 )- 



10 K.E. BASSLER, P.J. FORRESTER, AND N.E. FRANKEL 

The values of c\ and c-i follow by comparing (|3.24[) in the case p = 1. 

In the case [i = (|3.22[) and (I3.24|) reduce to first order recurrences. Their solutions imply 
the more explicit formulas (|3.17l) - (|3.18[) . □ 

3.3. The case a = 1. As revised below (|1.2[) . the case cr = 1 corresponds to the additive rank 1 
perturbation X ^ X + (/J,/(N + l))eief where X is a member of the (N + 1) x (N + 1) GUE. 
The correlation kernel for such random matrices has previously been expressed in the form [B] 

(3.26) K N+1 (x,y) = K% UE (x,y) + t^l e ^ 2 /2-y 2 /2 HN{x) f N+i ^. y) 

where 

f e -yz-z 2 /4 ^ z 

(3 ' 27) TN +^--=] C{o ^ } ^^ ) ^-- 

Here Cm. -2/4 is a simple closed contour encircling and — 2/i. Furthermore, we know p. 271) has 
the expansion [3 

(3.28) T N+1 (»,y)-(-l) {-j-yr - ^ {2jS)p+x ^.^ _ 1 _ p)x y 

To deduce (J3.26D from (|3.9p we first rearrange the latter to read 
K N+1 {x,y) = K^(x,y) 

+ e - ( ^ 2)/2 (tfW*) £ a N 0-p Nap 

\ =^-! NN-iyOCNPN-l ~ a N -iP N ) 

(3-29) +H N (x) ± hf ^-f~ aN - X \ M y)) ■ 

p—N 

Lemma 3.4. In the limit a — >• 1 

M v a N /3 p - p N a p { 1, p = N — 1 



yVjv-i cxatAv-i - fe«iv-i 10. P > N > 



Mn o<nI3n~i - Pnun- 



(2»y- N , p>n. 



Proof. These results follow by applying the recurrences Q3.22J) and (|3.24p to the numerators. 

□ 

Corollary 3.5. In the limit a — >• 1 

K N+1 (x,y) = 4 UE (^)+^ (lV)/2 %WE(2^ N ^T 

p=N yV P 

(3-30) = K^( X ,y) +e -(*W)/ 2 |^(>^ 2 - ^W^), 

where the second equality follows by taking <j — ¥ 1 in \S.15}) and substituting in the second equation 

of$m. 

Recalling p. 281) we see that (|3.30[) is identical to the known result (|3.26p . 
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3.6. The case \i = 0. Substituting (|3.17l) - (|3.19[) in (|3.8|) and (|3.9|) we obtain the specialization 
of the correlation kernel in the case [i — 0. 

Proposition 6. Suppose // = 0. For N even 

K N+1 (x,y) = K%Y$(x,y) 

gg-i(g) ^ 2 2 V(a 2 -l)Pg 2p+1 ( 2/ ) 

2^-2(^/2 _ 1)!(a 2 _ l)N /2-l ^Z._ i ^ p+1 

r », n , H N (x){N/2)\ ^ i (2p )! 2 

(3 - 31) + iV!(CT2 _ 1)J v/ 2p 2_, 2 A^^r ((T -^^m 

w/wte for N odd 
K N+1 {x,y) = K%™(x,y) 

, -tfW)l*( H N ^{x)((N-\)/2)\ ~ (2p)!(q 2 - l)Pff 2p (y) 

V(7V-l)!(a 2 -l)(^-D/ 2 p=( ^ i)/2 p!M p 

< 3 - 32 ) V-p-l^-lF-^ £ ^^(^-D-^xft,) 

vv II J \ I p=(JV-l)/2 p+ 

4. Soft edge scaling and the phase transition 

4.1. Secular equation. As revised in the Introduction, rank 1 perturbations can lead to a phase 
transition with respect to the location of the largest eigenvalue. In [3] we showed how for the 
additive rank 1 perturbation the corresponding secular equation associated with the eigenvalue 
problem could be used to predict the critical value s — l/v2JV (or equivalently fx = y/N/2) at 
which separation occurs, as well as the formula (jl.ip for the location of the largest eigenvalue 
in the separated phase. Here we will repeat those considerations in the case of the modified 
(N + 1) x (JV + 1) GUE matrices as specified by (|1.2p (the same reasoning applies to the GOE 
version of this construction). 

The secular equation determining the eigenvalues in this case is the condition (I2.5[) with {ftj} 
the eigenvalues of an N x N GUE matrix. Averaging over |v/v| 2 and cat it reduces to 



2 N i 



N 

A — a, 



3=1 



We seek the range of values of a 2 which permit this equation to be solved for |A| > y/2N, and the 
corresponding value of |A|. In view of the interlacing (I2.3|) this will correspond to the eigenvalues 
at the bottom and top of the spectrum separating from the GUE spectrum. We know the latter, 
to leading order, has density given by the Wigner semi-circle law 



\/9N / A 2 

( 4 - 2 ) ^)( A ) = VV 1 "^' |A|<2iV - 

But in general, for large N with {%} forming a continuum supported on / with density piu(y), 



3=1 



A — dj Jj\ — y 
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Substituting (|4.2j) in (|4.3I) the resulting integral can be evaluated (see e.g. [31 eq. (2.10)]). Thus 
/2JV »w 



%%*-*- f^)- — 



'-V2JV A-y * V V 2JV/' 
and so in the large AT limit the averaged secular equation (|3.17[) gives that for eigenvalue sepa- 
ration we must have 



(4.4) X -» = \ Y {l- ] Jl-— ), \\\>2N. 

Writing 



a 2 f A" \ 



(4.5) n = c^N/2 

and substituting A = y2N — the largest value of A which does not correspond to separation — 
we deduce that separation of the largest eigenvalue occurs for 

(4.6) a 2 + c> 2. 

Similarly, we deduce that separation of the smallest eigenvalue occurs for 

(4.7) a 2 - c> 2. 

And solving (|4.4[) under the conditions (I4.6[) and (|4.7p gives that the corresponding location of 
the separated eigenvalues will be at 

[W <x 4 + c 2 

1 ' ' V 2 (1 - a 2 /2)c± ct 2 ((c/2) 2 - (1 - a 2 )) 1 / 2 

with the choice of ± corresponding to (J4.6I) and (I4.7[) respectively. Note that when a = 1 this 
reclaims (jl.lj) . while for c = it reads 

IW 



(«) A = ± W2(a»-l)x/ a 

4.2. The case a = 1. The correlation kernel in this case is given by (13.301) or equivalently (13.26[) . 
We know from (jl.ll) that /i = y/N/2 is the critical value of /x in relation to separation of the 
largest eigenvalue from the Wigner semi-circle. A critical region occurs when the eigenvalues are 
scaled about the neighbourhood of the edge of the Wigner circle so that their spacing is 0(1), 
and the mean fi is scaled about the critical value yjN/2. Explicitly, in (|3.30p this is achieved by 
introducing the scaling variables X, Y, s according to 

(4.10) X = V2N +¥J ^, y = V2N +¥J ^, „ =y /Nj2(l + ^)- 
The task then is to compute the scaled N — > oo limit, which we know from [171 15] gives 

(4.11) y^—L—K N+1 (x t y) = K S °*(X,Y) + M(X) J e -l Y -»M(t)dt 

where K soit (X, Y) denotes the Airy kernel [S]. 

j^fv Y) = Ai(X)M'(Y) - M(Y)M'(X) 
1 ' ' X-Y 

As written the integral in (|4.1ip requires s > to be convergent, Noting that 

r-Y POO />00 

e -*(Y-t)Ai(t) dt = / e-< Y -^Ai(t)dt- / e^ s( - Y '^Ai(t)dt 



> J-oc JY 

/oo 
e -'<y-*)M{t)dt 



(4.13) 6 X 2H ,: {X) = 2 1 /4 n -Vi2 ( Ai(-«) + 0(n- 2 / 3 )) 
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gives this meaning for general s. The scaled correlation kernel (14.111) is precisely that which 
is found for the analogous scaling limit of the correlation kernel of the r = 1 spiked complex 
Wishart ensemble 0. 

4.3. The case /i = 0. The sums in Proposition[5]are very similar to those appearing in an earlier 
work [10J . There the main tool in the analysis of the soft edge limit of the sums is the uniform 
asymptotic expansion 

Proposition 7. Set 

a 2 = 2 - 2s/N 1/3 

and scale x and y as in \^. 10ty . Then fr4.11\ ) again results from the scaling limit of both the N 
even and N odd correlation kernels in Proposition 

Proof. Consider for definiteness the case N even. The terms in addition to K^^(x, y) of the 
correlation kernel can be written 

VV ; p=0 J *N-l+2p yV JV P =0 • /v Af+2p 

where 

_ / M N -i \i/* 2 2 r(N/2-l+p)l _ / Af N \ 1 / 2 (N/2y.(N + 2p)\ 

ap ' N ~ \Af N -i+2 P ) (N/2 -1)1 ' bp ' N ~W N+ap ) (N/2+p)\N\ ■ 

For large N and p ^C N we see from Stirling's formula that a Pt jq and b Pt N tend to 1. The key 
remaining step is to use (14.131) in the summands. Considering the first of these, for dchnitcncss, 
we substitute (|4~T^1) with n = N + 2p - 1 and -u = Y- 2p/N 1 / 3 . Writing 2pjN 1 ^ =iwe see 
that a Riemann integral approximation to the sum results, giving that the large N form of the 
latter is 

N l/4 roc 

(4.14) -^.J e- st Ai(Y-t)dt. 

Also, a literal application of (|4.13[) shows that the leading asymptotic form of the term outside 
the hrst sum is equal to 

(4.15) 2 1 / 4 iV- 1 / 12 AipO. 

The second term gives the same leading asymptotics. Thus we have to multiply twice (14.141) 
with flUTSJ) , and divide by 2 1 / 2 N 1 / 6 (recall the LHS of (ETTTj) in relation to this factor). The 
correction term to K soi \X, Y) on the RHS of (J4~TT|) results. □ 

4.4. General a 2 , fi > about eigenvalue separation. The correlation kernel (|3.29p con- 
sists of the GUE kernel and a correction term, the latter involving a summation over Hermite 
polynomials. The cases a 2 — 1 and /i = have tbe special feature that the coefficients exhibit 
simple functional forms (Corollary 13.51 and Proposition [5]). In contrast, for general a and \x the 
coefficients are given in terms of Hermite polynomials evaluated at a special point, and a certain 
Hilbcrt transformation of the same Hermite polynomial evaluated at this point. The key to cal- 
culating the large N, p, p <§; N form of the coefficients in this case relies on their characterization 
in terms of the solution of a difference equation, or more explicitly that of the scaled coefficients 

/a-ic\ ( -Miv-i \ 1//2 QW/3jv_i +p — PnQ.n~i+ p (i) 
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and 

(A 17 \ ( Mv \ 1/2 Pn&N-\+ p - ON0N-1+P _, ( 2 ) 

VA/"jv-i +p / /3at(57v-i — Av-icw 

where a p = Af p a p , j3 p = N p f3 p . 
Proposition 8. With \i replaced in favour of c according to C75l), and with 

c±(c 2 - 4(1 - er 2 )) 1/2 

(4.18) x± := ^^ ^ jj 

for N — > oo we have 

p p p p 

(4.19) y™ = X+X - X+X ~ , 7 < 2 > = ^J^-. 

y X- — x + ' x + — X- 

Proof. We know that both a p and (3 p satisfy the same recurrence Q3.22J) . But up to the scale 
factor (ATn-i+p) 1 2 , 7j> and j p are linear combinations of ajv-i+p and /3jv-i+p- Thus in 
(|3.22p . by replacing p by N — 1 + p, multiplying by through by (-/Vat-i+p) 1 ^ 2 and making the 

substitution (I4.5P we conclude that in the limit N — > oo both j p and y p satisfy the constant 
coefficient recurrence 

7-7+1 = c-lq + (o- 2 - 1)79-1 (9 = 1,2, ...)■ 

For 7p this is to be solved subject to the initial condition 70, 71 = 0, while for y p it is subject 
to the initial condition 70 = 0, 71 = 1. Solving the recurrence gives the stated result. □ 

Let us now write 

(4.20) c = £+T§7s, - 2 = (-) 2 Sl 



jyi/3' v > n 1 / 3 

with 
(4.21) (ct) 2 + c=2. 

According to (|4.6|) we are thus perturbing about the critical values for separation of the largest 
eigenvalue. Substituting in (J4.18I) shows 



(4.22) , W5 -l) + o(^) 



1 >\ (c-3)ai + 2s 2 

X ~ 2(c-2)iVV3 ' 



Restricting attention to the case /z > so that c > (recall the case c = as a special case 
covered above), and furthermore, from (|4.2ip . c < 2 (we exclude the case c = 2 since this would 
mean a 2 — which corresponds to a decoupled eigenvalue), we observe that \x+\ < 1. Recalling 
(|4.19p this tells us that for p large 

-v(i) _ x ~ ~(2) ^ (c- 1 )^- 

7 P c-2' 7p 8-2 ' 

Substituting in (|3.29p we therefore have that for N, p large, p <C N, 

(4.23) K N+1 (x,y)^K^(x,y)+e-^^^^^_^^.. 

•^N p=0 N p+N 

But this is precisely the setting of the sums in the proof of Proposition which we know give 
rise to the second term in (|4.1ip . 
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Proposition 9. Let c be related to fi by |^.5| ) ; scale c and a 2 according to l[4-2(fy with condition 
\4-21\) . Furthermore specify s by 

(4.24, * = ^|^ 

One has that the limiting correlation kernel is again given by (pTTTP- 

For c < we see from (|4.22[) that x+ < — 1. One is then faced with a formally divergent sum 

(4-25) £(_ 1)P | x+r ^WM, | X+ | > !. 

p=o N p+N 

Our method of analysis thus breaks down, whereby it was required that the limiting form of the 
coefficients be used inside of the summation. An alternative approach, with the sum rewritten 
— perhaps as an integral — is called for, although this is yet to be found. Nonetheless, with the 
case c = also giving (|4.11j) . we have no reason to expect anything different for c < 0, which 
would be the case if we could show that the term involving (J4.25I) vanished in the N — > oo limit. 
That the details of how this comes about are different may be due to the smallest (i.e. most 
negative) eigenvalue already being separated from the bulk of the spectrum in the case c < 0, 
but not in the case c > 0, an effect which must be encoded in the correction terms to the GUE 
kernel. 

This detail aside, our study then indicates that (|4.11[) is the universal correlation kernel for the 
scaling state about the separation of the largest eigenvalue due to a rank 1 type perturbation, as 
it is independent of the cause of the perturbation being with respect to the mean, the variance, 
or a combination of both. 

In Proposition [3] of Section 2 the eigenvalue PDF for the r-bordered GUE was given. To 
proceed as in Section 3 and give the explicit form of the correlation kernel would require inverting 
a 2r x 2r matrix (recall (J3.12I) '). This in turn puts a stop to us proceeding to study the scaling 
limit for general r. It is known in the special case a 1 = 1 from [171 |B]. 

It is natural to enquire into the analogous result for bordered GOE matrices. Very recently a 
method of analysis involving reduction to tridiagonal form [5] and stochastic differential equations 
has been introduced to study spiked real Wishart matrices [4]. A bordered GOE matrix can be 
reduced to tridiagonal form [12 . and one would expect this method to lead the conclusion that 
the eigenvalue separation effect for bordered GOE matrices shares the same universality class as 
that for spiked real Wishart matrices. 
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